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Abstract 

In this paper, we investigate dual 3-nets realizing the groups C3 x 
C3, C2 x C4, Alt4 and that can be embedded in a projective plane 
PG{2, K), where K is an algebraically closed field. We give a symbol- 
ically verifiable computational proof that every dual 3-net realizing 
the groups C3 x C3 and C2 x C4 is algebraic, namely, that its points 
lie on a plane cubic. Moreover, we present two computer programs 
qq ■ whose calculations show that the group Alt4 cannot be realized if the 

, characteristic of K is zero. 

cn 

: 1 Introduction 

In a projective plane a 3-net consists of three pairwise disjoint classes of lines 
such that every point incident with two lines from distinct classes is incident 
with exactly one line from each of the three classes. If one of the classes 
has finite size, say n, then the other two classes also have size n, called the 
order of the 3-net. In this paper we are considering 3-nets in a projective 
plane PG(2, K) over an algebraically closed field IK which are coordinatized 
by a group. Such a 3-net, with line classes A, B, C and coordinatizing group 
G = (G, •), is equivalently defined by a triple of bijective maps from G to 
(A,B,C), say 

a: G^A,f3: G ^ B, 7 : G^C 

such that a ■ b = c if and only if a (a), f3 (b), 7(c) are three concurrent lines 
in PG(2,K), for any a,b,c G G. If this is the case, the 3-net in PG(2,K) 
is said to realize the group G. Recently, finite 3-nets realizing a group in 
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the complex plane have been investigated in connection with complex line 
arrangements and resonance theory, see [1, 7, 9] and the references therein. 

Since key examples arise naturally in the dual plane of PG(2,K), it is 
convenient to work with the dual concept of a 3-net. Formally, a dual 3-net 
of order n in PG(2, K) consists of a triple (Ai, A 2 , A 3 ) with Ai, A 2 , A 3 pairwise 
disjoint point-sets of size n, called components, such that every line meeting 
two distinct components meets each component in precisely one point. A 
dual 3-net (Ai,A2, A3) realizing a group is algebraic if its points lie on a 
plane cubic. Moreover, we say that the dual 3-net (Ai,A 2 ,A 3 ) is of conic- 
line type ( triangular) if the components are contained in the union of a line 
and a nonsingular conic (in the union of three lines). 

In our computer-aided investigation, combinatorial methods are used to 
study finite 3-nets realizing the groups C 2 x C4, C 3 x C 3 , and Alt4. These 
results are fundamental for the complete classification of 3-nets embedded in 
a projective plane over a field, see [7]. Indeed, large groups could be dealt 
with theoretical results, but small groups having elements of order less than 
5 needed a more explicit computation. This is main motivation of this paper. 

We can summarize our results in the following theorem. 

Theorem 1.1. Let (A 1; A 2 , A 3 ) be a dual 3-net of order n which realizes a 
group G in the projective plane PG(2,K) defined over an algebraically closed 
field IK of characteristic p, where p = or p > 5. We also assume that n < p 
whenever p > 0. The following statements hold. 

(I) IfG = C 3 x C 3 or G = C 2 x C 4 , then (A u A 2 , A 3 ) is algebraic. 

(II) If p = 0, then the group Alt 4 cannot be realized. 

The proof is divided in three parts, see Sections 2,3, or 4, according as 
G = C 3 x C 3 , G = C 2 x C 4 , or G = Alt 4 . Our notation and terminology 
are standard, see [6]. In view of Theorem 1.1, K denotes an algebraically 
closed field of characteristic p where either p = or n < p where n denotes 
the order of the dual 3-net. 



2 



2 G = C 3 x C 3 



We denote by G = {0, . . . , 8} the elementary abelian group of order 9 given 
by the multiplication table 

012345678 

120453786 

201534867 

345678012 

453786120' 

534867201 

678012345 

786120453 

867201534 

Let H be the subgroup {0, 1, 2} of G. 

Let K be an algebraically closed field whose characteristic is either or 
more than 9. In this paper, all points are points of the projective plane over 
IK. We denote by oo a cubic root of unity in K. 

It is easy to see that realizations of the cyclic group of order 3 are pre- 
cisely the Pappus configurations. The point set V of a triangular dual 3-net 
realizing C3 consists of 9 points such that any line intersects V in 1 or 3 
points. In other words, V forms an AG(2,3), where AG(n, q) denotes the 
affme geometry over the finite field ¥ q . It is also well known that any AG (2, 3) 
embedded in PG(2, K) is a Hesse configuration, that is, the points are the 
inflection points of a nonsingular cubic curve. 

Lemma 2.1. Let 

A / = {0 1 ,1 1 ,2 1 ,0 2 ,1 2 ,2 2 ,0 3 ,1 3 ,2 3 } 

be a realization of H . Then there is a unique cyclic collineation a of order 
three mapping 

0i, li, 2i, 2 , 1 2 , 2 2 , 3 , 1 3 , 2 3 to 1 1; 2i, 0i, 1 2 , 2 2 , 2 , 2 3 , 3 , 1 3 , 

respectively, a is never central. The cubic curves containing A form a pencil. 
All these cubics are invariant under a. □ 
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Lemma 2.2. Let X be a set of nine points in a projective plane such that 
for all A,BeX, the line AB contains a third point of X . Then, X is either 
contained in a line, or form an AG{2,?>). □ 



In the sequel, we denote by A = {Oi, . . . , 81, O2, . . . , 82, O3, . . . , 83} a re- 
alization of G. We denote by A' the subset of A realizing the subgroup 
H = {0, 1, 2}. We will often use that the points of A can be re-indexed and 
the blocks {ii}, {32}, {^3} can be interchanged. 

Lemma 2.3. There is a line which intersects A in exactly two points. 

Proof. Assume that no line intersects A in exactly two points. As K* has 
no elementary abelian subgroup of order 9, A cannot be triangular or of 
conic-line type. Theorem 5.1 of [1] implies that none of the blocks {ii}, {j'2}, 
{k 3 } is contained in a line. Lemma 2.2 implies that these blocks must form 
an AG(2,3). Moreover, each line intersecting A in more than two points, 
intersect A in precisely three points. This means that with respect to the 
line intersections, A forms a Steiner triple system. As any three points of 
A generate a subsystem of order 9, A is in fact a Hall triple system, cf. [2, 
pages 496-499]. As | A| = 27, we obtain that A is an embedding of AG (3, 3) 
in a projective plane, which is not possible by [8] if char(K) 7^ 3. □ 

By re-indexing A, we can suppose that the line O1I1 intersects A in 
{Oi, li}, that is, Oi, li,2i are not collinear. Let a be the cyclic collineation 
of order three corresponding to the subnet A' realizing H = {0, 1,2}. We 
will choose our projective coordinate system such that the following hold: 

(1) d = (1, 0, 0),h = (0, 1, 0) and 2 l = (0, 0, 1). 

(2) F = (1, 1, 1) is a fixed point of a. 

(3) If the lines 0i0 2 , li I2, 2i2 2 are concurrent then F = (1,1,1) is their 
intersection. 

Notice that the lines contain no fixed point of a, hence (2) does not 
conflict with (1). Furthermore, if O1O2, I1I2, 2i2 2 are concurrent then their 
intersection is a fixed point of a. 

The collineation a has the matrix form 
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As O3 is not on the lines ljlj, O3 has coordinates of the form (a, 6,1) with 
a, 6 7^ 0. Then, we can compute the coordinates of the points I3 = (6, 1, a), 
2 3 = (1, a, b) and O2 = (6, 6a, a), I2 = (a, 6, 6a), 2 2 = (6a, a, 6). 

Let /3 be the cyclic collineation of order 3 corresponding to the subnet 
{0i, li, 2i, 3 2 , 4 2 , 5 2 , 3 3 , 4 3 , 5 3 }. The matrix of (5 has the form 

0m 
v 
1 

for some nonzero u, v G K. The point 3 2 has nonzero coordinates (x, 1). 
Then, we have 4 2 = (u,vx,y), 5 2 = (uy,uv,vx) and 

3 3 = (uy,vxy,vx), 4 3 = (xy,vx,y), 5 3 = (uvx,vuy,vxy). 

For all points ii, 2 G {3, . . . ,8}, there are three lines of the form j 2 /c 3 , 
j,k G {0, . . . , 5} such that i\ G j 2 /c 3 . The fact that the corresponding line 
triples are concurrent, can be expressed by the equations 

h = det(0 2 x 3 3 ,1 2 x 4 3 ,2 2 x 5 3 ) = 0, 

U = det(0 2 x 4 3 ,l 2 x 5 3 ,2 2 x 3 3 )= 0, 

h = det(0 2 x 5 3 ,1 2 x 3 3 ,2 2 x 4 3 ) = 0, 

U = det(3 2 x 3 ,4 2 x 1 3 ,5 2 x 2 3 ) = 0, 

h = det(3 2 x 1 3 ,4 2 x 2 3 ,5 2 x 3 ) = 0, 

/ 8 = det(3 2 x 2 3 ,4 2 x 3 ,5 2 x 1 3 ) = 0. 

The values a, 6, u, v, x, y determine A uniquely. The f^s (i G {3, . . . , 8}) are 
polynomial expressions of these values. In fact, we will look at a,b,u,v,x,y 
as indeterminates over IK and at the f^s as elements of K[o, 6, u, v, x, y]. The 
polynomials f7,fs,f9 have degree three in x,y, while for % = 4,5,6, the 
polynomials fi have the form /j = abvxyfi, where fi is in K[a, b, u, v, x, y\. 
The degree of f^, f 5 , / 6 in x, y is three. 

Generally speaking, we are looking for specializations such that the cor- 
responding configuration gives rise to a proper realization of G. 

Lemma 2.4. // any of the equations u — 1, v — 1, u — v holds then A is 
algebraic. 
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Proof. If u — v — 1 then a = (3. Let T be the cubic curve containing A' 
and 32- The equation of T can be computed explicitly, and one sees that if 
u — v — 1 then 33 G T. As by Lemma 2.1 T is invariant under ct = /3, we have 
42, 52, 4 3 , 53 G T, too. T cannot be completely reducible since then, some %\ 
would be collinear with some j 2 , k 2 . Suppose that r = £ U C with line £ and 
irreducible conic C. Then £, C are a-invariant and the lj's are in 3. If O2 G 3 
then O3, 13, 2 3 G £, and all 2,,'s are in 3 and all 3fc's are in £. As {Oi, li, 2i}, 
{0 2 ,l2,2 2 }, {0 3 ,1 3 ,2 3 }, {32,4 2 ,5 2 }, {3 3 ,4 3 ,5 3 } are all orbits of (a) and the 
lines O233, 1243, 2 2 53 are concurrent, we have that {3i, 4i, 5i}, {61, 7±, 81} are 
(a)-orbits contained in C . Continuing the process, we conclude that AcT 
(which is of course not possible). The same result is obtained if we start from 
j 2 G C or k 3 G C. 

Suppose now that V is irreducible. Denote by V* the set of nonsingular 
points. The (a)-orbits are all cosets of a subgroup H* of (T*,+) of order 
3. Then, simple arithmetic on T* yields that {3i,4i,5i}, {6i,7i,8i} are H* 
cosets of r*. Repeating this argument, we obtain A C T again. 

It remains to show that any of the equations u — 1, v — 1, u — v 
implies the other two. For that we observe the following equations of rational 
expressions: 
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In either case, the denominators at the left hand side cannot be zero as the 
corresponding lines are not concurrent. This proves that one equation implies 
another one, and two imply the third. This finishes the proof. □ 

The proof of the following lemma contains some elementary, but heavy 
computation. This computation can be formally verified by any computer 
algebra dealing with Groebner bases within a few seconds. 

Lemma 2.5. If a 3 = b 3 = 1 then v — 1. In particular, A is algebraic. 

Proof. We observe that a 3 = b 3 = 1 holds if and only if the lines O1O2, I1I2) 
2i2 2 and the lines Oil 2 , li2 2 , 2^2 are concurrent. (In other words, A' forms a 
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dual AG (2, 3).) Remember that in this case, our coordinate system is chosen 
such that F(l,l,l) is the fixed point 0i0 2 H Ul 2 H 2 X 2 2 . As d = (1,0,0) 
and O3 G 0\0 2 = 0\F, we have 6=1. By conjugaction in IK, we can assume 
a = oj w.l.o.g. 

Now, we can find polynomials Si,ti,pi,p 2 ,qi,q 2 , i £ {3,..., 8}, in the 
indeterminates a, 6, x, y, u, v with integer coefficients such that 

Sifi + pi(a — u) + ^2(6 — 1) = 18m; (t> — l)(t>:r — y 2 ){vx — uy 2 ), 
^2t i f i + q 1 (a-cj) + q 2 (b-l) = 18(v - l)(uy 3 - v 2 x 3 ). 

Assume v 7^ 1. Since det(0i, 3 2 , 4 2 ) = ?/ 2 - ra ^ 0, we have t>x — cuy 2 = 
wy 3 — f 2 x 3 = 0. This implies u = u 2 xy and v = uy 2 /x. Straightforward 
computation shows that the collineation 7 given by the matrix 

u 
cu 2 
1 

fixes 0i, li, 2i and maps the points O2, . . . , 52, O3, . . . , 53 to the points 

1 2 ,2 2 ,0 2 , 5 2 ,3 2 ,42, l3,2 3 ,0 3 , 53,3 3 ,4 3 , 

respectively. As 5 ■ 2 = 6 and 5 ■ 3 = 4 in G, we have 

7(3i) = 7(0 2 3 3 n 1 2 4 3 ) = 1 2 5 3 n 2 2 3 3 = 4 X . 

Similarly, 7(4^ = 5i and 7(5i) = 3i. Thus, 7 permutes the lines 3i0 3 , 
5i2 3 cyclically. As these lines intersect in 62, 62 is a fixed point of 7, which 
is not possible. □ 

We are now prepared to prove the main result. 

Theorem 2.6. A is algebraic. 

Proof. We can consider the /j's as polynomials in the indeterminates a, 6, 
u, v, x, y. Fix the values a,b,u,v and let Fi(X,Y) be the polynomials in 
two variables such that Fi(x,y) = fi(a,b,u,v,x,y). Define the linear series 
L generated by the Fj's. 

Recall that /3 is the collineation of order 2 mapping the points 0i, li, 2i, 
3 2 , 4 2 , 5 2 , 3 3 , 4 3 , 5 3 to U, 2i, 0i, 4 2 , 5 2 , 3 2 , 5 3 , 3 3 , 4 3 , respectively. From the 
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definition of the /j's one sees that the substitution X' — u/Y, Y' — vX/Y 
induces a linear automorphism of L of degree 3. We will denote this induced 
map by j3, as well. 

Define the polynomials 

E 1 = uvX + oo 2 uY 2 + uvX 2 Y, E 2 = vX 2 + u 2 uY + uY 2 X, 
E x = uvX + uuY 2 + co 2 vX 2 Y, E 2 = vX 2 + uuY + u 2 Y 2 X, 



and 



q 1 = uFl - F 8 , Q 2 = co 2 F A - F 5 , 

2 ; 



Qi = co 2 F 7 - F 8 , Q 2 = coF 4 - F t 
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of L. Then E 1: E 2 , Q\, Q 2 are eigenvectors of f3 with eigenvalue uuv/Y 3 and 
Ei, E 2 ,Q 1 ,Q 2 are eigenvectors of f3 with eigenvalue u 2 uv/Y 3 . We have the 
following resultant values: 

Re u e 2 (Y) = u 2 uvY(uv - Y 3 ) 2 ,R Bu e 2 (Y) = uuvY{uv - Y 3 ) 2 . 

This shows that the intersection of E\ — 0, E 2 — and the intersection of 
Ei — 0,Ei — consist of the points 0i, li, li (with multiplicity 0) and the 
fixed points of (5 (with multiplicity 2). In particular, E i: E 2 and Ei,E 2 are 
linearly independent. 

Straightforward calculation shows that 

Qi = GnEi + Gi 2 E 2 , Q 2 = GnEi + G 22 E 2 , 
Qi = GnEi + Gi 2 E 2 , Q 2 = GnEi + G 22 E 2 , 



where 



Gn = (u 2 b + ab 2 oj + a 2 )(u 2 + vu + u), 

Gi 2 = (b 2 co + tu 2 a + a 2 b){uv + tu 2 u + v to), 

G 2 i = u(u 2 b + ab 2 u + a 2 ){uj + u 2 v + u), 

G 22 = {b 2 u + u 2 a + a 2 b)(uu + u 2 v + uv), 

Gn = (ub + uj 2 ab 2 + a 2 )(u + u 2 v + u), 

Gi 2 = (au + a 2 b + u 2 b 2 )(uu + u 2 v + uv), 

G 2 i = tu 2 (cub + ijj 2 ab 2 + a 2 )(co 2 + v to + u), 

G 22 = (au + a 2 b + u 2 b 2 )(uv + u 2 u + vu>). 
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Assume that Q\,Q2 are linearly independent. Then Ei,E 2 G {Qi,Qq) < 
L. Therefore, Ti n • • • D T 4 is contained in the zero set of E\ — E% — 0, a 
contradiction. We can similarly show that Qi, Q2 must be linearly dependent. 
This implies 

= G11G22 — G12G21 

= (2 + u 2 )(b 2 u + u 2 a + ba 2 )(cuab 2 + u 2 b + a 2 ){u - v)(u - l)(v - 1), 

= G11G22 — G12G21 

= (2 + u)(au + ba 2 + u 2 b 2 )(ub + u 2 ab 2 + a 2 ){u - v)(u - l)(v - 1). 

The resultants of the polynomials (b 2 u + u 2 a + ba 2 )(uab 2 + u 2 b + a 2 ), (acu + 
ba 2 +tu 2 b 2 )(tub+uj 2 ab 2 +a 2 ) with respect to a, 6 are 96 7 (6 3 -l) 6 and 9a 7 (a 3 -l) 6 . 
Thus, A is algebraic by Lemmas 2.4 and 2.5. □ 

The Maple 13 program performing the computations of this section is at- 
tached in Appendix A. We use Buchberger's algorithm in order to explicitely 
construct the polynomials of Lemma 2.5. Thus, any computer algebra which 
can do symbolical calculation with rational polynomials can be used to verify 
the results. This convinces us about the correctness of our computation. 



3 G = C 2 x C 4 

The main ingredient of the proof is Lame's Theorem [7, Proposition 2.3]. 
A classical Lame configuration consists of two triples of distinct lines in 
PG(2, K), say £±, £2, £3 and r 1? r 2 , r 3 , such that no line from one triple passes 
through the common point of two lines from the other triple. For 1 < j, k < 3, 
let Rjk denote the common point of the lines £j and r^. There are nine such 
common points, and they are called the points of the Lame configuration. 

Proposition 3.1. Lame's Theorem. If eight points from a Lame configura- 
tion lie on a plane cubic then the ninth also does. 
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The group C 2 x C 4 can be given by the multiplication table 

1 2 3 4 5 6 7 8 

2 1 4 3 6 5 8 7 

3 4 1 2 7 8 5 6 

4 3 2 1 8 7 6 5 

5 6 7 8 2 1 4 3' 

6 5 8 7 1 2 3 4 

7 8 5 6 4 3 2 1 

8 7 6 5 3 4 1 2 

The triple {ii,j 2 , ks} is collinear if and only if i * j = k. 

The following 6-tuples of collinear points form a Lame configuration: 

C/x = {U, 1 2 , 1 3 }, {3!, 5 2 , 7 3 }, {6!, 7 2 , 3 3 } + {li, 7 2 , 7 3 }, {3!, 1 2 , 3 3 }, {6!, 5 2 , 1 3 }, 
C/ 2 = {li, 3 2 , 3 3 }, {3!, 7 2 , 5 3 }, {6i, 5 2 , 1 3 } + {lx, 5 2 , 5 3 }, {3 1; 3 2 , 1 3 }, {6!, 7 2 , 3 3 }, 
U 3 = {l u 1 2 , 1 3 }, {3x, 7 2 , 5 3 }, {8 1} 5 2 , 3 3 } + {lx, 5 2 , 5 3 }, {3 1; 1 2 , 3 3 }, {8 1} 7 2 , 1 3 }, 
U 4 = {h, 3 2 , 3 3 }, {3!, 5 2 , 7 3 }, {8i, 7 2 , 1 3 } + {li, 7 2 , 7 3 }, {3 1; 3 2 , 1 3 }, {8i, 5 2 , 3 3 }. 

Let C be a cubic curve through the points 

li) I2, 13, 3i, 7 2 , 7 3 , 61, 7 2 , 3 3 . 

Then \C fl U\\, \C fl U 2 \ > 8, hence, C passes through the nineth points 3 2 
and 5 2 . It follows that \C D U 3 \, \C (1 U 4 \ > 8. Thus, C contains 

U ± U U 2 U U 3 U U A = {li, 3i, 61, 81, 1 2 , 3 2 , 5 2 , 7 2 , 1 3 , 3 3 , 5 3 , 7 3 }. 

It is straightforward to check that any of the following Lame configurations 
intersects C in at least 8 points: 

{lx, 5 2 , 5 3 }, {3x, 1 2 , 3 3 }, {5x, 3 2 , 7 3 } + {l u 3 2 , 3 3 }, {3x, 5 2 , 7 3 }, {5 1; 1 2 , 5 3 }, 
{li, 1 2 , 1 3 }, {3i, 5 2 , 7 3 }, {7i, 3 2 , 5 3 } + {lx, 5 2 , 5 3 }, {3i, 3 2 , 1 3 }, {7 U 1 2 , 7 3 }, 
{lx, 5 2 , 5 3 }, {61, 7 2 , 3 3 }, {81, 2 2 , 7 3 } + {lx, 7 2 , 7 3 }, {61, 2 2 , 5 3 }, {8 1? 5 2 , 3 3 }, 
{lx, 5 2 , 5 3 }, {61, 4 2 , 7 3 }, {81, 7 2 , 1 3 } + {li, 7 2 , 7 3 }, {6 1? 5 2 , 1 3 }, {8 1? 4 2 , 5 3 }, 
{lx, 1 2 , 1 3 }, {61, 7 2 , 3 3 }, {81, 3 2 , 6 3 } + {lx, 3 2 , 3 3 }, {61, 1 2 , 6 3 }, {81, 7 2 , 1 3 }, 
{lx, 1 2 , 1 3 }, {61, 3 2 , 8 3 }, {81, 5 2 , 3 3 } + {lx, 3 2 , 3 3 }, {6 1; 5 2 , 1 3 }, {8 1} h, 8 3 }. 
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Hence, C contains the further points 5i, 7i, 2 2 , 4 2 , 6 3 , 8 3 . Finally, we consider 
the Lame configurations 

{lx, 1 2 , 1 3 }, {2 X , 5 2 , 63}, {61, 2 2 , 5 3 } + {li, 5 2 , 5 3 }, {2 X , 2 2 , 1 3 }, {6!, 1 2 , 63}, 
{3i, 1 2 , 3 3 }, {4 1; 7 2 , 63}, {6!, 2 2 , 5 3 } + {3i, 7 2 , 5 3 }, {4 X , 2 2 , 3 3 }, {6 1? 1 2 , 63}, 
{lx, 5 2 , 5 3 }, {7!, 6 2 , 3 3 }, {8 1? 3 2 , 63} + {li, 6 2 , 63}, {7 X , 3 2 , 5 3 }, {8 1? 5 2 , 3 3 }, 
{li, 8 2 , 83}, {5!, 3 2 , 7 3 }, {6!, 7 2 , 3 3 } + {li, 7 2 , 7 3 }, {5!, 8 2 , 3 3 }, {6!, 3 2 , 83}, 
{lx, 1 2 , h}, {7i, 7 2 , 2 3 }, {8!, 2 2 , 7 3 } + {lx, 2 2 , 2 3 }, {7 1? 1 2 , 7 3 }, {8 1? 7 2 , 1 3 }, 
{3x, 2 2 , 4 3 }, {5x, 1 2 , 5 3 }, {6x, 7 2 , 3 3 } + {3x, 1 2 , 3 3 }, {5 1; 7 2 , 4 3 }, {6 1? 2 2 , 5 3 }. 

As before, one sees that any of them has at least 8 points in common with 
C, thus, C passes through all the points of the embedding of C 2 x C4. 



4 G = Alt 4 (p = 0) 



Let the group Alt (4) be given on the underlying set {1, . . . , 12} by the Cayley 
table 
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projective plane form a realization of Alt (4), if for all k = 1, . . . , 12, 
iii 3ii are collinear if and only if i * j — k. 

Proposition 4.1. Alt(4) cannot be realized on the complex projective plane. 

Proof. We see that {1, 2, 3, 4} is an elementary Abelian normal subgroup and 
{1, 5, 9} is a subgroup. Without a loss of generality, we can assume that 

li = [1,0,0], 1 2 = [0,1,0], i 3 = [1,-1,0], 

2 1 = [0,1,1], 2 2 = [1,0,1], 2 3 = [0,0,1]. 
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As 3i, 1 2 , 3 3 are collinear, we can take 3i and 3 3 in the form 3i = [a, b, c] and 
3 3 = [a, b + 1, c]. This enables us to compute the remaining points 

3i = [a,b,c], 3 2 = [a — 1, 1 + b, c], 3 3 = [a, 6 + 1, c], 

4i = [a- l,l + 6,c- 1], 4 2 = [a,6,c - 1], 4 3 = [a - 1, b, c - 1]. 

These points indeed form a realization of C 2 x C 2 . Similarly, we choose 
5i and 9i as generic points [di,d 2 ,l], [^4,^5,1]. Then 5 3 ,9 3 have the form 
[di, d 3 , 1], [d 4 , d 6 , 1], respectively and computation yields 

5 2 = [^4 + d 5 - d 3 , d 3 , 1], 9 2 = [di + d 2 - d 6 , d 6 , 1]. 

Define the points 

5i = 9 2 l 3 n 1 2 5 3 , 5 2 = 9il 3 n li5 3 , 5 3 = 5il 2 n li5 2 , 

61 = 9 2 4 3 n 2 2 5 3 , 6 2 = 9 1 2 3 n 4 1 5 3 , 6 3 = 5i2 2 n 4^2, 

7x = 9 2 2 3 n 3 2 5 3 , 7 2 = 9x3 3 n 2^3, 7 3 = 5x3 2 n 2^2, 

81 = 9 2 3 3 n 4 2 5 3 , 8 2 = 9i4 3 n 3i5 3 , 8 3 = 5i4 2 n 3i5 2 , 

9i = 5 2 l 3 n 1 2 9 3 , 9 2 = 5il 3 n li9 3 , 9 3 = 9il 2 n li9 2 , 

10i = 5 2 3 3 n 2 2 9 3 , 10 2 = 5i2 3 n3i9 3 , 10 3 = 9i2 2 n 3i9 2 , 

Hi = 5 2 4 3 n 3 2 9 3 , 11 2 = 5 1 3 3 n 4 1 9 3 , 113 = 9x32 4192, 

12i = 5 2 2 3 n 4 2 9 3 , 12 2 = 5x4 3 n 2x9s, 12 3 = 9^ n 2 ± 9 2 . 

Let us denote by dijk the determinant of the 3x3 matrix with rows i±, j 2 , k 3 . 
We define the sets X = {d ijk | k = i * j} and 

Y = {(1,2,6), (9,1,10), (1,10,8), (5,9,2), (9,5,2), (1,10,2), 

(1.3.11) , (1,4,12), (1,5,7), (9,1,12), (9,5,4), (5,9,4), 

(1.5.6) , (1,6,11), (1,9,12), (1,12,4), (1,5,8), (1,7,3), 

(1.7.12) , (5,9,3), (9,1,11), (1,12,7), (1,9,11), (5,1,6), 

(5.1.7) , (1,6,2), (1,3,7), (1,8,4), (1,11,6), (5,1,8), 
(9,5,3), (1,8,10), (1,11,3), (1,4,8), (1,9,10), (1,2,10)}. 

For all (i,j,k) e Y, i • j ^ k, thus for any proper complex realization of 
Alt (4), we can substitute complex numbers in the variables a, b, c, di, ■ • ■ ,d$ 
such that all polynomials in X are zero and all for all k) G Y, d^k 7^ 0. 
Put g = IL( i>jtk)eY dijk and define 

X> = {f/gcd(f,g)\feX}. 
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It is still true that Alt (4) has a realization if and only if one can substitue 
complex numbers in a, b, c, d±, ■ ■ • , d$ such that all polynomials in X' vanish. 
Groebner basis computation shows that the ideal generated by X' contains 
the polynomials d\ — d^, d2 — d%, implying 5i = 9i, a contradiction. □ 

The Groebner basis computation of this section is too heavy for most of 
computer programs of this type. We found two programs which is able to 
compute the Groebner basis: the F4 algorithm [4] in the computer algebra 
system Maple 13 and the modStd library [5] of Singular [3]. These pro- 
grams do not store the cofactors of the Groebner bases, hence one cannot 
verify the result symbolically. However, two completely different implemen- 
tations deliver the same result, thus, we can trust this computation as well. 

The Maple 13 implementation is attached in Appendix B and the Sin- 
gular implementation is attached in Appendix C. The computations take 
less that 3 minutes, and less that 3 hours, respectively. 

A Maple code for the case G = C3 x C3 

This appendix contains 

######################################################## 

# Maple 13 program for computing with dual 3-nets 

# G = C3 x C3 

######################################################## 

# Preparation 

with(LinearAlgebra) : 

isect : =proc (a,b, c ,d) 

evala(CrossProduct(CrossProduct(a,b) .CrossProduct (c ,d) ) ) : 
end proc : 

idet : =proc (a, b , c , d , e , f ) 

evala(Determinant(<CrossProduct(a,b) I CrossProduct (c ,d) I CrossProduct (e ,f )>) ) : 
end proc : 

alias (omega=RootOf (X"2+X+l) ) ; 

######################################################## 

# Part 1: Constructing the points, the transformations and the equations. 

# One defines the transformations, the base points and the equations which 

# correspond to certain collinearities . The unknowns $a,b,u,v,x,y$ are seen as 

# fixed elements of the base field. 

# In the program, Px_i denotes the point $x_i$ of the dual $3$-net. 

alpha : =«0, 1,0> I <0,0, 1> I <1,0,0>>; 
beta:=«0,v,0>|<0,0,l>|<u,0,0»; 
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P0_1:=<1,0,0>: 
P3_2:=<x,y, 1>; 
P0_3:=<a,b,l>; 



Pl_l 
P4_2 
Pl_3 



=<0,1,0>: P2_1:=<0,0,1>: 
=beta.P3_2; P5_2:=beta.P4_2; 
= (alpha"(-D) .P0_3; P2_3 : = (alpha" (-1) ). Pl_3 ; 
# We turn Lambda_3 in the opposite direction! 



P0_2 : =isect (P0_1 ,P0_3 , Pl_l , Pl_3) ; 
Pl_2:=isect(P0_l,Pl_3,Pl_l,P2_3); 
P2_2 : =isect (P0_1 ,P2_3 ,P1_1 ,P0_3) ; 
P3_3 : =isect (P0_1 ,P3_2 ,P1_1 ,P5_2) ; 
P4_3 : =isect (P0_1 ,P4_2 ,P1_1 ,P3_2) ; 
P5_3 : =isect (P0_1 ,P5_2 ,P1_1 , P4_2) ; 



f:=[ 0,0, 



idet (P0. 


_2,P3. 


_3,P1_ 


_2,P4. 


_3,P2_ 


_2,P5_ 


_3)/(a*b*v*x*y) , 


# 


P3_ 


.1 


idet(P0. 


.2,P4. 


_3,P1. 


.2,P5. 


.3,P2. 


.2,P3. 


_3)/(a*b*v*x*y) , 


# 


P4_ 


.1 


idet(P0. 


.2,P5. 


_3,P1. 


.2,P3. 


.3,P2. 


.2,P4. 


_3)/(a*b*v*x*y) , 


# 


P5_ 


.1 


idet (P3_ 


.2,P0. 


.3,P4. 


.2, PI. 


.3,P5. 


_2,P2. 


.3), 


# 


P6. 


.1 


idet (P3. 


.2, PI. 


.3,P4. 


.2,P2. 


.3,P5. 


_2,P0_ 


.3), 


# 


P7. 


.1 


idet(P3. 


_2,P2. 


_3,P4_ 


_2,P0_ 


_3,P5. 


.2, PI. 


.3) 


# 


P8_ 


.1 



] : 



f :=f actor (f) : 

######################################################## 

# Part 2 (Lemma 2.4): Any of u=l, v=l, u=v implies the other two equations. 

factor (subs (u=l , f [6] /idet (P0_2 , P4_3 , P3_2 , P0_3 , P4_2 , Pl_3) ) ) ; 
factor (subs (v=l , f [6] /idet (P0_2 , P5_3 , Pl_2 , P3_3 , P3_2 , P0_3) ) ) ; 
factor (subs (u=v, f [6] /idet (P0_2 , P3_3 , Pl_2 , P4_3 , P3_2 , P0_3) ) ) ; 

######################################################## 

# Part 3 (Lemma 2.5): If a~3=b~3=l then v=l. 

# We can assume a=omega and b=l. 

gb : =Groebner [Basis] ( [op(f ) ,a-omega,b-l] , plex(u, v,x,y ,a,b) , output=extended) : 
f actor (gb[l] [3]), factor (gb[l] [5]); 

# We construct the cof actors explicitely: 
s:=evala(18*gb[2] [3] [1. .8]) : 
t:=evala(18*gb[2] [5] [1. .8]) : 
q:=evala(18*[gb[2] [5] [9] ,gb[2] [5] [10]]): 
p:=evala(18*[gb[2] [3] [9] ,gb[2] [3] [10]]): 

factor (add(f [i]*s[i] ,i=3. . 8)+p [1] * (a-omega)+p [2] * (b-1) ) ; 
factor (add (f [i] *t [i] ,i=3. . 8)+q[l] * (a-omega)+q[2] * (b-1) ) ; 

# This shows that all cof actors have coefficients in Z [omega] : 
seq(denom(f actor (_u) ) ,_u in [op(s) ,op(t) ,op(p) ,op(q)] ) ; 

######################################################## 

# Part 4: Computation with the beta-invariant polynomials. 

# From now on, we consider $a,b,u,v$ as fixed elements of the base field 

# and $X,Y$ as indetereminates . 

# We define the action of $\beta$ on the polynomial ring in two variables. 
F:=map(_x->subs({x=X,y=Y},_x) ,f ) : map(_x->degree(_x,{X,Y}) ,F) ; 
betaonpoly : =proc (U) return factor (subs ({X=u*Y/ (v*X) ,Y=u/X} ,U) ) : end proc : 
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# This shows that the nontrivial solutions of F [3] =F [4] =F [5] =0 

# and F [6] =F [7] =F [8] =0 are $\beta$-invariant : 

map(_x->f actor(_x) , 

[bet aonpoly (F [3] ) /F [5] , bet aonpoly (F [4] ) /F [3] , betaonpoly (F [5] ) /F [4] ] 

); 

map(_x->f actor(_x) , 

[betaonpoly (F [6] )/F [7] .betaonpoly (F [7] ) /F [8] .betaonpoly (F [8] ) /F [6] ] 

); 

# We define the E[i]'s, barE[i]'s, Q[i]'s and barQ[i]'s 

# and show that their curves are $\beta$-invariant : 

E:=[u*v*X+omega"2*u*Y"2+omega*v*X~2*Y, v*X"2+omega"2*u*Y+omega*Y"2*X] : 
barE:=[u*v*X+omega*u*Y"2+omega~2*v*X~2*Y,v*X~2+omega*u*Y+omega"2*Y"2*X] : 

Q : = [omega * F [6] -F [7] , omega" 2 * F [3] -F [4] ] : 
barQ : = [omega" 2 * F [6] -F [7] , omega * F [3] -F [4] ] : 

seq(f actor(betaonpoly(_u)/_u) , _u in [op(E) ,op(Q)] ) ; 
seq(f actor(betaonpoly(_u)/_u) , _u in [op(barE) ,op(barQ)] ) ; 

# We define the G[i,j]'s and barG[i,j]'s. 

# We show that they are indeed coefficients of Q[i]'s and barQ[i]'s. 

G [2, 1] : =omega* (omega"2*b+a*b"2*omega+a"2) * (omega+omega"2*v+u) ; 
G[2,2] :=(b"2*omega+omega"2*a+a"2*b)*(omega*u+omega"2*v+u*v) ; 
G[l,l] :=(omega"2*b+a*b"2*omega+a"2)*(omega"2+v*omega+u) ; 
G[l,2] : =(b"2*omega+omega"2*a+a"2*b) * (u*v+omega"2*u+v*omega) ; 
barG[l,l] : =(omega*b+omega"2*a*b"2+a"2) * (omega+omega"2*v+u) ; 
barG [1,2] : =(omega*a+a"2*b+omega"2*b"2) * (omega*u+omega"2*v+u*v) ; 
barG [2,1]: =omega"2* (omega*b+omega"2*a*b"2+a"2) * (omega"2+v*omega+u) ; 
barG [2 , 2] : =(omega*a+a"2*b+omega"2*b"2) * (u*v+omega"2*u+v*omega) ; 

map(_x->evalb(f actor(_x) ) , [ 
Q[1]=G[1,1]*E[1]+G[1,2]*E[2] , 
q [2] =G [2 , 1] *E [1] +G [2 , 2] *E [2] , 
barQ [1] =barG [1,1] *barE [1] +barG [1,2] *barE [2] , 
barQ [2] =barG [2 , 1] *barE [1] +barG [2 , 2] *barE [2] 

]); 

# We compute the factors of the determinants of the G[i,j]'s and barG[i,j]'s: 

map(_x->evalb(f actor(_x) ) , [ 
G [1 , 1] *G [2 , 2] -G [1 , 2] *G [2,1] = 

(2+omega"2)*(b"2*omega+omega"2*a+b*a"2)*(omega*a*b"2+omega"2*b+a"2)*(u-v)*(u-l)*(v-l) , 
barG [1,1] *barG [2 , 2] -barG [1,2] *barG [2,1] = 

(2+omega)*(a*omega+b*a"2+omega"2*b"2)*(omega*b+omega"2*a*b"2+a"2)*(u-v)*(u-l)*(v-l) 

]); 

ra:=resultant( 

(-b+a*b"2*omega-omega*b+a"2) * (-b"2*omega+a+a*omega-b*a"2) , 
(-omega*b+a*b"2+a*b"2*omega-a"2)* (a*omega+b*a"2-b"2-b~2*omega) , 
a): 

rb:=resultant( 

(-b+a*b"2*omega-omega*b+a"2) * (-b"2*omega+a+a*omega-b*a"2) , 
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(-omega*b+a*b~2+a*b~2*omega-a~2)* (a*omega+b*a~2-b~2-b~2*omega) , 
b): 

factor(ra/(b~3-l)~6) ; 
factor(rb/(a-3-l)~6) ; 



B Maple code for the case G = Alt4 

This appendix contains the implementation of the computations of Section 
4, using the F4 algorithm [4] in the computer algebra system Maple 13. 
This program does not store the cofactors of the Groebner bases, hence one 
cannot verify the result symbolically. The computation takes less than 3 
minutes. 

######################################################## 

# Maple 13 program for computing with dual 3-nets 

# G = Alt (4) 

######################################################## 

# Preparation 

with (Linear Algebra) ; 
isect :=proc(a,b, c,d) 

evala(CrossProduct(CrossProduct(a,b) .CrossProduct (c ,d) ) ) : 
end proc: 

ct : =Matrix( 

[ [ 1,2,3,4,5,6,7,8,9,10,11,12 ], 

[ 2,1,4,3,7,8,5,6,12,11,10,9 ], 

[ 3,4,1,2,8,7,6,5,10,9,12,11 ], 

[ 4,3,2,1,6,5,8,7,11,12,9,10 ], 

[ 5,6,7,8,9,10,11,12,1,2,3,4 ] , 

[ 6,5,8,7,11,12,9,10,4,3,2,1 ] , 

[ 7,8,5,6,12,11,10,9,2,1,4,3 ] , 

[ 8,7,6,5,10,9,12,11,3,4,1,2 ] , 

[ 9,10,11,12,1,2,3,4,5,6,7,8 ] , 

[ 10,9,12,11,3,4,1,2,8,7,6,5 ] , 

[ 11,12,9,10,4,3,2,1,6,5,8,7 ] , 

[ 12,11,10,9,2,1,4,3,7,8,5,6 ] 

]); 

d:=array(l . .6) ; 

######################################################## 

# Part 1: We define the points of the dual 3-net 

# using a,b,c,d[l] , . . .d[6] as indeterminates . 

P:=[ [<1,0,0>,<0,1,0>,<1,-1,0>] , 
[<0,1,1>,<1,0,1>,<0,0,1>] , 
[<a,b,c>,0,<a,l+b,c>] , 
[0,0,0] , 

[<d[l] ,d[2] ,l>,0,<d[l] ,d[3] ,1>] , 
[0,0,0] , 
[0,0,0] , 
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[0,0,0] , 

[<d[4] ,d[5] ,l>,0,<d[4] ,d[6] ,1>] , 
[0,0,0] , 
[0,0,0] , 
[0,0,0] ]; 

# As P[4,l], P[l,2], P[4,3] are coll, we may assume wlog that 

# P[4,l]=<a,b,c> and P [4,3] =<a, 1+b, c> . 

# Similar argument for P[5,3] and P[9,3], using the fact that 

# these points cannot have last coordinate 0. 

P[3,2] :=evala(isect(P[3,l] ,P[1,3] ,P[1,1] ,P[3,3])/c); 
P[4,2] :=evala(isect(P[3,l] ,P[2,3] ,P[2,1] ,P[3,3])/a) ; 

P[4,l] :=evala(isect(P[3,2] ,P[2,3] ,P[2,2] ,P [3 , 3] )/ (1+b) ) ; 
P[4,3] :=evala(isect(P[l,l] ,P[4,2] ,P[2,1] ,P [3 , 2] )/ (1+b-c) ) ; 

P[5,2] :=isect(P[l,l] ,P [5 , 3] ,P [9 , 1] ,P [1 ,3] ) ; 
P[9,2] :=isect(P[l,l] ,P [9 ,3] ,P [5 , 1] ,P [1 ,3] ) ; 

######################################################## 

P[5,l] :=isect(P[9,2] ,P [1 ,3] ,P [1 , 2] ,P [5 ,3] ) : 
P[5,2] :=isect(P[9,l] ,P [1 ,3] ,P [1 , 1] ,P [5 ,3] ) : 
P[5,3] :=isect(P[5,l] ,P[1,2] ,P[1,1] ,P[5,2]): 

P[6,l] :=isect(P[9,2] ,P[4,3] ,P[2,2] ,P[5,3]) : 
P[6,2] :=isect(P[9,l] ,P[2,3] ,P[4,1] ,P[5,3]) : 
P[6,3] :=isect(P[5,l] ,P[2,2] ,P[4,1] ,P[5,2]): 

P[7,l] :=isect(P[9,2] ,P [2 , 3] ,P [3 , 2] ,P [5 ,3] ) : 
P[7,2] :=isect(P[9,l] ,P [3 ,3] ,P [2 , 1] ,P [5 ,3] ) : 
P[7,3] :=isect(P[5,l] ,P[3,2] ,P[2,1] ,P[5,2]): 

P[8,l] :=isect(P[9,2] ,P[3,3] ,P[4,2] ,P[5,3]) : 
P[8,2] :=isect(P[9,l] ,P [4 ,3] ,P [3 , 1] ,P [5 ,3] ) : 
P[8,3] :=isect(P[5,l] ,P[4,2] ,P[3,1] ,P[5,2]): 

P[9,l] :=isect(P[5,2] ,P [1 , 3] ,P [1 , 2] ,P [9 ,3] ) : 
P[9,2] :=isect(P[5,l] ,P [1 , 3] ,P [1 , 1] ,P [9 ,3] ) : 
P[9,3] :=isect(P[9,l] ,P[1,2] ,P[1,1] ,P[9,2]): 

P[10,l] :=isect(P[5,2] ,P [3,3] ,P [2, 2] ,P [9 , 3] ) : 
P[10,2] :=isect(P[5,l] ,P [2 ,3] ,P [3, 1] ,P [9 ,3] ) : 
P[10,3] :=isect(P[9,l] ,P[2,2] ,P[3,1] ,P[9,2]) : 

P[ll,l] :=isect(P[5,2] ,P[4,3] ,P[3,2] ,P[9,3]) : 
P[ll,2] :=isect(P[5,l] ,P[3,3] ,P[4,1] ,P[9,3]) : 
P[ll,3] :=isect(P[9,l] ,P [3,2] ,P [4, 1] ,P [9 ,2] ) : 

P[12,l] :=isect(P[5,2] ,P [2 ,3] ,P [4, 2] ,P [9 , 3] ) : 
P[12,2] :=isect(P[5,l] ,P [4,3] ,P [2, 1] ,P[9,3] ) : 
P[12,3] :=isect(P[9,l] ,P [4,2] ,P [2, 1] ,P [9 ,2] ) : 



######################################################## 

# Part 2: We construct the polynomial identities. 
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eqs : = [] : 

for i from 1 to 12 do 
for j from 1 to 12 do 

aa: determinant (<P[i,l] I P [j ,2] I P [ct [i , j] ,3] >) : 
eqs := [op (eqs) ,aa] : 
end do 
end do : 

######################################################## 

# Part 3: We filter out the nonzero factors of the equations. 

nepos : = [ 

[1, 2, 6], [9, 1, 10], [1, 10, 8], [5, 9, 2], 
[9, 5, 2], [1, 10, 2], [1, 3, 11], [1, 4, 12], 
[1, 5, 7], [9, 1, 12], [9, 5, 4], [5, 9, 4], 
[1, 5, 6], [1, 6, 11], [1, 9, 12], [1, 12, 4], 
[1, 5, 8], [1, 7, 3], [1, 7, 12], [5, 9, 3], 
[9, 1, 11], [1, 12, 7], [1, 9, 11], [5, 1, 6], 
[5, 1, 7], [1, 6, 2], [1, 3, 7], [1, 8, 4], 
[1, 11, 6], [5, 1, 8], [9, 5, 3], [1, 8, 10], 
[1, 11, 3], [1, 4, 8],[1, 9, 10], [1, 2, 10] 
]; 

noneqs:=map(_x->Determinant(<P[_x[l] ,1] |P[_x[2] ,2] |P[_x[3] ,3]>) , nepos) : 

noneqs:=mul(x,x in noneqs) : 

eqs : =select (x->x<>0 , eqs) : nops (eqs) ; 

eqs_reduced : =map(x->f actor (x/gcd(x , noneqs) ) , eqs) : 

map (degree , eqs) -map (degree , eqs_reduced) ; 

######################################################## 

# Part 4: We compute the Groebner basis of the corresponding ideal. 

# This Groebner basis shows that d[l]=d[4], d[2]=d[5], d[3]=d[6]. 

# The computation takes less that 3 minutes using the F4 algorithm. 

gb:=Groebner [Basis] (eqs_reduced, tdeg(a,b, c ,d [1] ,d[2] ,d[3] ,d[4] ,d[5] ,d[6])); 



C Singular code for the case G = Alt4 

This appendix contains the implementation of the computations of Section 
4, using the modStd library [5] of Singular [3]. This program does not 
store the cofactors of the Groebner bases, hence one cannot verify the result 
symbolically. The computation takes less than 3 hours. 

//////////////////////////////////////////////////////// 

// Singular 3.1 program for computing with dual 3-nets 
// G = Alt (4) 

//////////////////////////////////////////////////////// 
// Preparation 

LIB "modstd.lib" ; 
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intmat ct [12] [12] = 

I, 2,3,4,5,6,7,8,9,10,11,12 , 
2,1,4,3,7,8,5,6,12,11,10,9 , 
3,4,1,2,8,7,6,5,10,9,12,11 , 
4,3,2,1,6,5,8,7,11,12,9,10 , 

5.6.7.8.9.10.11.12.1.2.3.4 , 

6.5.8.7.11.12.9.10.4.3.2.1 , 
7,8,5,6,12,11,10,9,2,1,4,3 , 

8.7.6.5.10.9.12.11.3.4.1.2 , 
9,10,11,12,1,2,3,4,5,6,7,8 , 

10.9.12.11.3.4.1.2.8.7.6.5 , 

II, 12,9,10,4,3,2,1,6,5,8,7 , 
12,11,10,9,2,1,4,3,7,8,5,6; 

ring r=0, (a,b,c,d(l. .6)) ,dp; 

//////////////////////////////////////////////////////// 
// Part 1: We define the points of the dual 3-net 
// using a,b,c,d[l] , . . .d[6] as indeterminates . 

// As P[4,l], P[l,2], P[4,3] are coll, we may assume wlog that 
// P[4,l]=<a,b,c> and P [4,3] =<a, 1+b, c> . 

// Similar argument for P[5,3] and P[9,3], using the fact that 
// these points cannot have last coordinate 0. 



list pt_data= 



1,0,0, 


0,1,0, 


1,-1,0, 


// 


1 


0,1,1, 


1,0,1, 


0,0,1, 


// 


2 


a,b,c, 


0,0,0, 


a, l+b,c, 


// 


3 


0,0,0, 


0,0,0, 


0,0,0, 


// 


4 


d(l),d(2),l, 


0,0,0, 


d(l),d(3),l, 


// 


5 


0,0,0, 


0,0,0, 


0,0,0, 


// 


6 


0,0,0, 


0,0,0, 


0,0,0, 


// 


7 


0,0,0, 


0,0,0, 


0,0,0, 


// 


8 


d(4) ,d(5) ,1, 


0,0,0, 


d(4) ,d(6) ,1, 


// 


9 


0,0,0, 


0,0,0, 


0,0,0, 


// 


10 


0,0,0, 


0,0,0, 


0,0,0, 


// 


11 


0,0,0, 


0,0,0, 


0,0,0; 


// 


12 



//////////////////////////////////////////////////////// 

// Procedures for manipulating the points of the dual 3-net: 

proc rpt(int x, int i) 
{ 

return (list (pt_data[(x-l)*9+(i-l)*3+l. . (x-1) *9+(i-l) *3+3] ) ) ; 

} 

proc setpoint(int x, int i, list u) 
{ 

pt_data[(x-l)*9+(i-l)*3+l]=u[l] ; 
pt _data [ (x- 1 ) *9+ ( i- 1 ) *3+2] =u [2] ; 
pt .data [ (x- 1 ) *9+ ( i- 1 ) *3+3] =u [3] ; 

} 

proc divpoint(int x, int i, poly p) 
{ 
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pt_data[(x-l)*9+(i-l)*3+l]=pt_data[(x-l)*9+(i-l)*3+l] /p; 
pt_data[(x-l)*9+(i-l)*3+2]=pt_data[(x-l)*9+(i-l)*3+2] /p; 
pt_data[(x-l)*9+(i-l)*3+3]=pt_data[(x-l)*9+(i-l)*3+3] /p; 

} 

proc crossprodQist u,list v) 
{ 

return(list (u [2] *v [3] -u [3] *v [2] , -u [1] *v [3] +u [3] *v [1] , u [1] *v [2] -u [2] *v [1] ) ) ; 

} 

proc isect(u,v,w,z) 
{ 

return(crossprod(crossprod(u,v) , crossprod(w.z) ) ) : 

} 

proc detpoints(int x, int y, int z) 
{ 

matrix m [3] [3] = 

pt_data[(x-l)*9+0*3+l. . (x-1) *9+0*3+3] , 

pt_data[(y-l)*9+l*3+l. . (y-1) *9+l*3+3] , 

pt_data[(z-l)*9+2*3+l. . (z-1) *9+2*3+3] ; 
return(det(m)) ; 

} 

//////////////////////////////////////////////////////// 
// Part 2: We set the remaining points of the dual 3-net. 

setpoint(3,2,isect(rpt(3,l) ,rpt(l,3) ,rpt(l,l) ,rpt(3,3))) ; 
divpoint(3,2,c) ; print (rpt (3,2) ) ; 

setpoint(4,2,isect(rpt(3,l) ,rpt(2,3) ,rpt(2,l) ,rpt(3,3))) ; 
di vpoint (4 , 2 , a) ; print (rpt (4,2)); 

setpoint(4,l,isect(rpt(3,2) , rpt (2, 3) , rpt (2, 2) , rpt (3, 3))) ; 
di vpoint (4 , 1 , 1+b) ; print (rpt (4,1)); 

setpoint(4,3,isect(rpt(l,l) ,rpt(4,2) ,rpt(2,l) ,rpt(3,2))) ; 
divpoint(4,3, 1+b-c) ; print (rpt (4, 3)) ; 

//////////////////////////////////////////////////////// 
setpoint(5,2,isect(rpt(l,l) ,rpt(5,3) ,rpt(9,l) ,rpt(l,3))) ; 
setpoint(9,2,isect(rpt(l,l) ,rpt(9,3) ,rpt(5,l) ,rpt(l,3))) ; 

//////////////////////////////////////////////////////// 
setpoint(5,l,isect(rpt(9,2) ,rpt(l,3) ,rpt(l,2) ,rpt(5,3))) ; 
setpoint(5,2,isect(rpt(9,l) ,rpt(l,3) ,rpt(l,l) ,rpt(5,3))) ; 
setpoint (5,3, isect (rpt (5, 1) , rpt (1 ,2) ,rpt (1 , 1) ,rpt (5,2)) ) ; 

setpoint(6,l,isect(rpt(9,2) ,rpt(4,3) ,rpt(2,2) ,rpt(5,3))) ; 
setpoint (6, 2, isect (rpt (9,1) ,rpt(2,3) ,rpt(4,l) ,rpt(5,3))) ; 
setpoint (6, 3, isect (rpt (5,1) ,rpt(2,2) ,rpt(4,l) ,rpt(5,2))) ; 

setpoint (7,1, isect (rpt (9, 2) , rpt (2, 3) , rpt (3, 2) , rpt (5, 3))) ; 
setpoint (7, 2, isect (rpt (9,1) ,rpt(3,3) ,rpt(2,l) ,rpt(5,3))) ; 
setpoint (7, 3, isect (rpt (5,1) ,rpt(3,2) ,rpt(2,l) ,rpt(5,2))) ; 

setpoint (8,1, isect (rpt (9, 2) , rpt (3, 3) , rpt (4, 2) , rpt (5, 3))) ; 
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setpoint(8,2,isect(rpt(9,l) ,rpt(4,3) ,rpt(3,l) ,rpt(5,3))) ; 
setpoint(8,3,isect(rpt(5,l) ,rpt(4,2) ,rpt(3,l) ,rpt(5,2))) ; 

setpoint(9,l,isect(rpt(5,2) ,rpt(l,3) ,rpt(l,2) ,rpt(9,3))) ; 
setpoint(9,2,isect(rpt(5,l) ,rpt(l,3) ,rpt(l,l) ,rpt(9,3))) ; 
setpoint(9,3,isect(rpt(9,l) ,rpt(l,2) ,rpt(l,l) ,rpt(9,2))) ; 

setpoint(10,l,isect(rpt(5,2) ,rpt(3,3) ,rpt(2,2) ,rpt(9,3))) ; 
setpoint(10,2,isect(rpt(5, 1) ,rpt(2,3) ,rpt(3,l) ,rpt(9,3))) ; 
setpoint(10,3,isect(rpt(9, 1) ,rpt(2,2) ,rpt(3,l) ,rpt(9,2))) ; 

setpoint(ll,l,isect(rpt(5,2) ,rpt(4,3) ,rpt(3,2) ,rpt(9,3))) ; 
setpoint(ll,2,isect(rpt(5, 1) ,rpt(3,3) ,rpt(4,l) ,rpt(9,3))) ; 
setpoint(ll,3,isect(rpt(9, 1) ,rpt(3,2) ,rpt(4,l) ,rpt(9,2))) ; 

setpoint(12,l,isect(rpt(5,2) ,rpt(2,3) ,rpt(4,2) ,rpt(9,3))) ; 
setpoint(12,2,isect(rpt(5, 1) ,rpt(4,3) ,rpt(2,l) ,rpt(9,3))) ; 
setpoint(12,3,isect(rpt(9, 1) ,rpt(4,2) ,rpt(2,l) ,rpt(9,2))) ; 



//////////////////////////////////////////////////////// 
// Part 3: We define the nonzero polynomials. 

intmat nz_pos [36] [3] = 



1,2,6, 


9,1,10, 


1,10,8, 


5,9,2, 


9,5,2, 


1,10,2, 


1,3,11, 


1,4,12, 


1,5,7, 


9,1,12, 


9,5,4, 


5,9,4, 


1,5,6, 


1,6,11, 


1,9,12, 


1,12,4, 


1,5,8, 


1,7,3, 


1,7,12, 


5,9,3, 


9,1,11, 


1,12,7, 


1,9,11, 


5,1,6, 


5,1,7, 


1,6,2, 


1,3,7, 


1,8,4, 


1,11,6, 


5,1,8, 


9,5,3, 


1,8,10, 


1,11,3, 


1,4,8, 


1,9,10, 


1,2,10; 



list nz; 

for (int i=l; i<=36; i++) 
{ 

nz=insert (nz , detpoints (nz_pos [i , 1] , nz_pos [i , 2] , nz_pos [i , 3] ) ) ; 

} 

proc nonzero_reduction(poly p) 
{ 

for (int i=l; i<=36; i++) 
{ 

p=p/gcd(p,nz[i] ) ; 

} 

return (p) ; 

} 

//////////////////////////////////////////////////////// 
// Part 4: We construct the polynomial identities. 

poly p; 
ideal 1=0; 

for (int i=l; i<=12; i++) 
{ 

for (int j=l; j<=12; j++) 
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{ 

p=detpoints(i , j ,ct [i , j] ) ; 

if (p!=0) { I=nonzero_reduction(p) ,1; } 

} 

} 

//////////////////////////////////////////////////////// 
// Part 5: We compute the Groebner basis of the corresponding ideal. 
// This Groebner basis shows that d[l]=d[4], d[2]=d[5], d [3] =d [6] . 
// The computation takes less that 3 hours using the modStd method. 

ideal J=modStd(I) ; 

J; 
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